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1.0  SUMMARY 


The  objective  of  this  study  was  to  assess  the  feasibility  of  develop¬ 
ing  a  single  crystal  pressure  sensor  using  Stoneley  waves.  A  theoretical 
study  of  Stoneley  wave  propagation  for  the  case  of  a  single  material  but  dif¬ 
ferent  orientations  was  carried  out  on  cubic  materials  as  well  as  non-cubic 
but  piezoelectric  quartz  and  lithium  niobate,  LiNbO-j.  The  results  of  these 
studies  show  that  strongly  hound  Stoneley  waves  propagate  along  the  major  axes 
of  cubic  materials.  Similar  analysis  carried  out  on  quartz  and  LiNbO-j  show 
that  only  the  YZ  cut  of  UNb03  supports  Stoneley  waves.  These  waves  are  loosely 
hound  and  have  an  electromechanical  coupling  coefficient  of  1  »  10"^.  Our 
conclusions  are  that  a  Stoneley  wave  sensor  is  feasible  on  YZ  LifJbC^  but  not 
on  the  major  cuts  of  quartz.  We  recommend  continued  studies  to  optimize 
Stoneley  wave  search  techniques,  studies  of  more  materials  for  Stoneley  wave 
devices,  and  experimental  device  fabrication  to  confirm  these  results. 
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2.0  BACKGROUND 

High  sensitivity,  underwater  sensors  are  needed  for  acoustic  surveil¬ 
lance  at  great  depths,  typically  1000-3000  ft.  Hydrostatic  pressures  at  these 
depths  are  extremely  large  reaching  values  as  high  as  1500  psi.  Existing 
ceramic  transducers  using  the  piezoelectric  effect  experience  depoling,  aging, 
and  other  undesirable  behavior  under  large  hydrostatic  pressures.  To  avoid 
depoling  and  to  enhance  their  sensitivity,  such  sensors  are  used  with  a  high 
dc  bias  voltage.  This  leads  to  complicated  and  cumbersome  electronic  arrange¬ 
ments.  Many  of  these  ceramic  materials  are  sensitive  to  water,  thus  requiring 
costly  protective  coatings.  All  of  these  problems  compound  an  already  severe 
problem  of  instability  or  aging  of  the  ceramic  material  itself.  There  is  at 
present  a  need  for  underwater  sensors  which  are  capable  of  high  sensitivity 
under  large  hydrostatic  pressure,  low  in  cost,  rugged,  insoluble  in  seawater, 
and  require  simple  electronics. 

The  objective  of  this  program  was  to  assess  the  feasibility  of  de¬ 
veloping  single  crystal  pressure  sensors  for  towed  array  applications  under 
large  hydrostatic  pressures  using  a  Stoneley  wave  structure  as  shown  in 
Fig.  1.  Our  goal  was  to  conduct  a  theoretical  and  if  possible  experimental 
study  of  Stoneley  wave  propagation  and  sensors.  Hence  a  major  part  of  this 
program  was  to  develop  a  computer  program  to  study  the  wave  propagation  of 
Stoneley  waves  in  anisotropic  piezoelectric  materials. 

The  operation  of  the  Stoneley  wave  sensor  shown  is  similar  to  a 
recently  developed^  surface  acoustic  wave  hydrophone  sensor.  In  these 
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structures  a  high  Q  resonating  electrode  pattern  is  used.  Flexure  and/or 
compression  of  the  crystal  causes  a  frequency  shift  from  the  quiescent 
resonant  frequency.  If  the  crystal  is  used  as  the  frequency  control  element 
of  a  highly  stable  oscillator,  frequency  modulated  (FM)  sidebands  are  created 
when  underwater  sound  impinges  on  the  sensing  crystal.  Using  surface  acoustic 
wave  resonators,  sensitivities  as  high  as  -ISO  dB  re  1  V/u  Pa  have  been  demon¬ 
strated.  1 

Surface  acoustic  wave  sensing  crystals  cannot  be  used  under  large 
hydrostatic  pressures  because  only  one  side  of  the  crystal  can  be  exposed  to 
the  surrounding  seawater.  However  the  Stoneley  wave  structure  does  not  have 
this  limitation.  Because  the  Stoneley  waves  are  bound  to  the  interface  be¬ 
tween  two  dissimilar  crystals  and  do  not  interact  with  the  outward  surfaces, 
the  resonant  electrode  pattern  is  isolated  *rom  these  surfaces.  This  fact 
would  allow  hydrostatic  pressure  to  be  equalized  by  applying  water  to  each 
side  of  the  structure. 

To  date  Stonely  wave  have  not  been  used  in  resonant  crystal  struc- 

p 

tures.  In  general  two  dissimilar  materials  are  required11  to  sustain  Stoneley 
waves  at  their  common  boundary.  Not  all  materials,  however,  form  the  proper 
boundary  conditions.  In  fact,  only  a  very  few,  a  point  made  by  a  studyJ  of 
900  materials  in  which  only  30  combinations  resulted  in  Stoneley  waves. 

However,  as  pointed  out  by  lim  and  Musgrave,^  a  single  crystalline 
material  can  be  used  for  both  top  and  bottom  materials,  provided  different 
orientations  are  used.  The  boundary  in  this  case  if  formed  due  to  the 
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anisotropic  nature  of  the  interface.  Lim  and  Musgrave  demonstrated  this  using 
single  crystal  copper,  a  simple  non-piezoelectric  cubic  material.  In  this 
report  their  work  has  been  extended  to  more  complex  anisotropic  crystals  which 
are  also  piezoelectric. 
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3.0  THEORETICAL  FORMULATION 


3. 1  Stonele.y  Wave  Propagation 

Consider  an  interface  wave  propagating  along  the  X-axis  of  two 
anisotropic  media  as  shown  in  Fig.  2.  The  X-axis  is  normal  to  the  interface 
itself.  Any  wave  propagating  in  this  system  must  satisfy  the  wave  equation. 


(1) 


where  T^  is  the  stress  tensor,  Uj  is  the  mechanical  displacement  and  p  is  the 
density.  The  electrical  displacement  current  density,  D,  for  a  dielectric 
material  must  also  satisfy  the  divergence  relationship. 


3D. 

_ 7_ 

3Xi 


(2) 


For  piezoelectric  materials,  the  equations  of  state  are  as  follows: 


Tij  =  ci  jkl  Skl  '  ek i j  Ek 


(3) 
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(b) 


Fig.  2  (a)  Surface  acoustic  wave  boundary  conditions; 

(b)  Stoneiey  wave  boundary  conditions. 
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-ikl  °kl 


+  e 


ik  k 


(4) 


where  the  strain  tensor,  ,  is  related  to  components  of  displacement  by 


S 


(5) 


and  the  electric  field,  E,  is  expressed  as  the  negative  gradient  of  a  scalar 
potential , 


(6) 


The  elastic  constants  of  the  material  are  expressed  by  the  tensor,  c^kl;  the 
piezoelectric  properties  by  the  tensor,  e^j;  and  the  dielectric  properties  by 
the  tensor,  eik.  In  this  formulation  the  piezoelectric  elements  are  evaluated 
at  constant  electric  field  and  dielectric  properties  are  evaluated  at  constant 
strain.^ 

Substituting  Eqs.  (3)-(6)  into  (1)  and  (2)  yields  the  following 

equations 
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3U 
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3U 
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i) 

H 


d2p 

ekij  JTTX^  = 


d  p 


2 
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(7) 

(8) 


A  Stoneley  wave  propagating  along  the  interface  is  assumed  to  have  the  form. 


i  kbX-j  ik(x.-vt) 

Uj  -  e  J  e  1  (9) 


where  Uj  are  particle  displacement  vectors  along  Xj-axes;  wavenumber 
k  =  2  Tt/X;  b  is  the  decay  constant  with  depth;  and  v  is  the  phase  velocity. 
The  scalar  electric  potential  is  assumed  to  follow  a  similar  relationship, 


i  k  bX ,,  ik(x.-vt) 
$  =  a4  e  J  e  1 


(10) 


Substituting  these  propagating  wave  forms  into  the  above  equations  of 
motion  ((7), (8))  yields  the  secular  matrix  equation, 
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where  =  055^  +  2ci5b  + 

r22  =  c44b2  +  2c46b  +  c66 

r33  3  c33b2  +  2c35b  +  c55 

r 12  3  c45b2  +  <c14  +  c56^b  +  c16 

r13  =  c35b2  +  (c13  +  c55^b  +  c15  (12) 

r23  3  c34b2  +  (c36  +  c45^b  +  c56 

r14  3  e35b2  +  (e15  +  e3l)b  +  ell 

r24  3  e34b2  +  ^e14  +  ®36^b  +  e16 

r34  =  e33b2  +  (e13  +  e35^b  +  e15 

and  the  electric  term  is  given  by 

r44  =  -(c33b2  +  2ei3b  +  en)  (13) 
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The  secular  Eq.  (11)  is  now  applied  to  the  top  and  bottom  substrate 
regions  and  the  matrix  of  coefficients  evaluated.  This  yields  a  matrix  of 
coefficients  in  the  unknown  decay  constant  b  and  the  wave  velocity  v.  For  a 
given  wave  velocity  v,  the  matrix  of  coefficients  is  expanded  into  a  poly- 
nominal  in  b  which  in  the  completely  general  case  is  an  8th  order  polynomial. 
The  solutions  to  the  polynomial  are  8  complex  roots  or  decay  constants,  bn. 
The  roots  occur  normally  in  conjugate  pairs  and  only  those  roots  which  lead  to 
decay  with  distance  away  from  the  interface  in  each  region  are  considered. 

This  procedure  yields  a  total  of  8  roots,  4  from  each  region  of  the  interface. 
Substituting  these  roots  back  into  the  secular  equation  then  yields  the  rela¬ 
tive  vectors  for  the  particle  displacements,  ajn^,  and  the  potential,  a|n). 

From  the  above  analysis  it  is  clear  that  8  partial  solutions  are 
required  to  completely  satisfy  any  and  all  boundary  conditions.  Thus  the 
total  solution  to  the  layered  problem  is  taken  as  the  sum, 


5j  =  £  Cn  aj  e 
T,8 


(n)  JkbnX3  jMXrvt) 


(14) 


■E 

T  ,8 


,  ( n)  JkbnX3  Jk<Xr¥t> 


cn  a4  e 


(15) 


11 

C3015A/bw 


Rockwell  international 


ERC41040. 1FR 

where  T,B  indicates  summation  over  root  numbers,  n,  for  the  top  and  bottom 
substrates.  Each  solution  is  scaled  by  the  appropriate  complex  constant,  Cn, 
corresponding  to  the  appropriate  complex  decay  constant,  bn. 

The  complete  solution  is  obtained  by  satisfying  eight  boundary 
conditions.  These  boundary  conditions  are  as  follows: 

a.  Transverse  Mechanical 


1.  Continuity  of  transverse  displacement  at  the  interface, 


2.  Continuity  of  transverse  stress  at  the  interface, 


b.  Electrical  Continuity 


1.  Continuity  of  the  normal  component  of  displacement  current  at 
the  interface, 
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2.  Continuity  of  the  scalar  potential  at  the  interface. 


c.  Sagittal  Plane  Mechanical  Displacements 


1.  Continuity  of  the  longitudinal  component  of  mechanical 
displacement  at  the  interface. 


2.  Continuity  of  the  vertical  component  of  mechanical  displacement 
at  the  interface, 


3.  Continuity  of  the  shear  component  of  the  stress  at  the 
interface. 


4.  Continuity  of  the  compressional  component  of  the  stress  at  the 
interface. 
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If  the  Interface  is  conducting  then  the  conditions  on  the  electric 
displacement  current  and  potential  are  changed  to  reflect  a  net  zero  potential 
at  the  interface  for  both  the  top  and  bottom  regions. 

Substituting  the  summation  solutions  (14),  (15)  into  the  above 
boundary  conditions  results  in  a  set  of  8  homogeneous  equations  in  the  unknown 
constants,  Cn-  The  matrix  coefficients  are  shown  in  detail  in  Table  I.  The 
determinant  is  commonly  referred  to  as  the  boundary  condition  determinant.  In 
order  to  have  a  non-trivial  solution  the  determinant  must  equal  zero.  A  com¬ 
puter  search  routine  is  used  to  determine  when  and  if  such  a  velocity  exists. 
After  the  Stoneley  wave  velocity  has  been  found  (not  always  possible),  the 
matrix  of  boundary  condition  coefficients  is  used  to  determine  the  relative 
scaling  vectors,  Cn,  and  the  solution  is  complete. 

3. 2  Program  Software  Development 

A  completely  general  set  of  computer  programs  were  written  to  obtain 
Stoneley  wave  solutions  in  terms  of  the  top  and  bottom  material  piezo-elastic 
and  dielectric  material  constants.  These  constants  were  obtained  from  well 
known  studies6  of  material  properties.  Since  the  material  constants  are 
always  given  with  respect  to  the  major  axes  of  the  crystal,  a  subprogram  was 
used  to  rotate  the  constants  to  the  appropriate  off-axes  directions  to  be 
studied.  After  these  rotations  a  search  was  performed  to  determine  if  a 
Stoneley  wave  solution  existed. 
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Table  I 


Stone! ey  Wave  Boundary  Condition  Determinant 


Bottom 

Substrate  Columns 
m  (Four  Values) 
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m  (Four  Values) 
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The  Stone! ey  wave  search  routine  was  based  upon  the  value  of  the 
determinant  of  the  boundary  condition  equations.  This  determinant  was 
evaluated  as  a  function  of  velocity,  a  solution  being  indicated  wherever  the 
value  of  the  determinant  was  zero.  For  each  velocity  the  secular  equation, 

Eq.  (11),  was  determined  and  the  partial  solution  vectors  and  associated  decay 
constants  evaluated  using  a  convergent  algorithm^  for  solving  the  eighth  order 
polynomi nal . 

Whenever  a  Stoneley  wave  solution  was  found  to  exist,  an  iteration  on 
velocity  using  a  gradient  search  routine  was  performed.  This  technique  was 
able  to  determine  the  correct  velocity  to  nine  significant  digits.  After 
determining  the  velocity  to  the  desired  accuracy,  the  displacement  and  poten¬ 
tial  fields,  Eq.  (14)- (1 5) ,  were  then  plotted  as  a  function  of  X3  away  from 
the  interface. 


► 
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4.0  RESULTS 


The  Stone! ey  wave  formulation  and  resulting  computer  programs  des¬ 
cribed  in  Section  3.0  were  used  to  study  propagation  in  a  number  of  materials 
considered  to  be  readily  available  should  an  appropriate  Stoneley  wave  solu¬ 
tion  be  found.  Initially  our  studies  concentrated  on  verifying  the  existence 
of  Stoneley  waves  in  simple  non-piezoelectric  cubic  materials  and  determining 
search  techniques.  After  these  studies  more  complicated  crystal  systens  such 
as  quartz  and  LiNbOj  were  examined.  A  considerable  amount  of  computer  "data" 
was  accumulated  and  in  the  interest  of  keeping  this  report  readable  only  a 
representative  amount  will  be  presented. 

4. 1  Cubic  Crystals 

As  reported  by  lim  and  Musgrave  a  single  crystal  material  with  dif¬ 
ferent  orientations  for  top  and  bottom  materials  can  support  a  Stoneley 
wave.  To  be  described  in  this  section  are  some  prototypic  results  using  the 
slightly  modified  set  of  elastic  constants  for  copper  shown  in  Table  II.  The 
specific  modification  involved  only  C44  which  was  changed  such  that  the 
anisotropy  ratio 


Y 


2C 


44 


C1I  '  C12 


1.52 
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Table  II 

Modified  Copper  Elastic  Constants 
(*  1011  Nt/m3) 


This  modification  reflected  a  ratio  which  more  closely  paralleled  the  ratio 
found  in  quartz  and  LiNbO-j  crystals.  The  general  character  of  the  Stoneley 
wave  solutions  was  not  substantially  changed  by  the  modification. 

Rotated  elastic  matrices  for  Stoneley  waves  were  generated  using 
contra-rotated  elastic  constants  for  the  top  and  bottom  media  as  shown  in 
Fig.  3.  In  this  case  rotation  was  in  the  plane  of  the  interface.  Shown  in 
Table  III  are  the  rotated  elastic  matrices  for  +5°  copper. 


Table  III 

Rotated  Copper  Elastic  Constants 
(*  1011  Nt/m2) 


1.714 

1.235 

1.239 

• 

• 

1 0. 0206 

1.235 

1.714 

1.239 

• 

• 

+0. 0206 

1.239 

1.239 

1.71 

• 

• 

• 

• 

• 

• 

0.  356 

• 

« 

• 

• 

• 

• 

0.356 

• 

iO. 0206 

?0. 0206 

• 

• 

• 

0. 3525 
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Fig.  3  Rotation  of  material  matrices  and  vectors  for  Stoneley 
wave  formulation. 
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After  rotating  the  material  constants,  a  Stoneley  wave  search  was 
conducted.  For  simple  cubic  materials  such  as  copper,  Stoneley  waves  were 
always  found  to  exist  over  a  broad  range  of  angular  rotations.  Solution 
details  for  four  different  angles  of  rotation  are  tabulated  in  Tables  IV- 
VII.  In  each  table  the  complex  secular  equation  roots  (decay  constants), 
partial  solution  displacement  vectors,  end  complex  solution  weighting 
constants  are  tabulated.  The  total  solution  displacements  on  either  side  of 
the  interface  are  plotted  in  Figs.  4-7. 

Examination  of  Figs.  4-7  as  well  as  Table  IV-VII,  reveals  character¬ 
istics  which  were  found  to  be  typical  of  all  Stoneley  wave  solutions: 

1.  The  Stoneley  wave  velocity  was  always  numerically  greater  than 
the  equivalent  surface  wave  velocity  but  lower  than  the  slowest 
shear  wave  velocity. 

2.  Small  angular  material  rotations  always  produce  solutions  with 
small  decay  constants,  i.e.,  the  wave  is  not  well  bound  to  the 
interface. 

3.  Large  angular  material  rotations  always  produce  well  bound 
Stoneley  waves  with  predominantly  transverse  shear  displace¬ 
ments,  eventually  degenerating  into  a  pure  shear  displacement 
wave. 
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Fig.  6  Stoneley  wave  displacements  for  copper  ±10°. 
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Table  IV 

Stoneley  Wave  Solution  for  ±2°  Copper 
Velocity  =  1959.152  m/s 


Roots 

bN 

Displacement  Vectors 

ai 

Weighting  Factors 

CN 

00 

f  4 

o 

1 

o 

• 

o 

a1  =  1.0  +  jO.O 
<*2  =  64.  736  ♦  jO.  0 
a3  =  0.0  -  j0.476 

0.0  -  jO.  008 

0.0  +  jO. 184 

=  1.0  ♦  jO.O 
a.£  =  64.  736  +  jO.  0 
<*3  «  0.0  +  J6.476 

0.0  +  jO. 008 

0.403  ♦  jO.  007 

aj  =  1.0  ♦  jO.O 
a2  3  -0.113  -  jO. 003 
a3  *  2.192  -  jO. 032 

1.0  +  jO.O 

r-* 

,  o 

1  ° 

k 

!  - 

CO 

o 

o 

aj  =  1.0  +  jO.O 
a2  =  0. 118  -  jO. 003 
a3  =  -2.192  *  jO. 032 

1 

00 

o 

o 

O  I 

•n 

» 

o 

*—«  1 
1 

P  ! 
-C* 

o 

00 

I 

t-,.  1 

o 

o 

o 

t  . 

aj  3  l.  J  *  jO.O 
a 2  3  -0. 118  +  jO. 003 
a3  =  2.192  -  jO.  032 

1.0  +  jO.  008 

-0.403  +  jO. 007 

ax  =  1.0  +  jO.O 
a 2  3  0. 1 18  +  jO. 003 
a3  =  2. 192  ♦  jO. 032 

1.0  +  jO.O 
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Table  V 

Stoneley  Wave  Solution  for  ±5°  Copper 
Velocity  =  1956.870  m/s 


N 

Roots 

bN 

Displacement  Vectors 

ai 

Weighting  Factors 

CN 

1 

0.0  +  jO. 170 

=  1.0  +  jO. 0 
a2  -  -32.  747  +  jO. 0 
a3  =  0. 0  +  j7. 514 

-0.001  +  jO. 044 

2 

0.0  -  jO. 170 

=  1.0  +  jO.O 
<*2  ~  32.  747  +  jO.  0 
a-j  =  0.0  -  j 7 .  514 

O 

O 

o 
►— + 

1 

o 

o 

4* 

- - — 

3 

O 

-t* 

o 

4^ 

1 

O 

o 

-c* 

CO 

_ 

aj  =  1.0  +  jO.  0 

c*2  =  -0.  297  -  jO.  055 
a3  =  -2.186  -  jO. 212 

O 

O 

+ 

o 

f—* 

1 

1 

0.404 

aj  3  1.0  +  jO.O 
a2  3  0.297  -  jO. 055 
a3  =  -2.136  +  jO. 212 

-0.999  -  jO. 044 

5 

-0.404  -  jO. 043 

dj  =  1.0  +  jO.O 
a2  *  0. 297  +  jO. 055 
a3  =  2.186  -  jO. 212 

0.999+  jO.  044 

6 

-0.404  +  jO. 043 

a1  =  1.0  +  jO.O 
a2  *  0. 297  +  jO. 055 
a3  «  2.186  +  jO. 212 

1.0  +  jO.O 
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Table  VI 

Stoneley  Wave  Solution  for  tlO°  Copper 
Velocity  =  1936.134  m/s 


N 

Roots 

bN 

Displacement  Vectors 

ai 

Weighting  Factors 

CN 

1 

0.0  +  jO. 140 

<*!  =  1.0  +■  jO.O 
<%2  ~  62. 177  +  jO. 0 
a3  =  0.0  +  j 1 6. 912 

-0.001  -  jO. 051 

2 

0.0  -  jO. 140 

aj  =  1.0  +  jO.O 
<*2  =  62.  177  +  jO.  0 
a3  =  0.0  -  j 16. 912 

0.001  -  jO.  051 

3 

0.410  -  jO. 164 

al  •  i-o  ♦  JO.O 
<*2  s  -0.408  -  jO.  354 
a3  *  1.9333  -  jO. 697 

-1.0  +  jO.O 

1 

0.401  +  jO. 164 

*  1.0  +  jO.O 
a 2  =  0.408  -  jO. 354 
a3  =  1.933  +  jO. 697 

-0.999  -  jO. 051 

5 

-0.401  +  jO. 164 

<*1  =  1.0  +  jO.O 
=  0.408  +  j0.354 
a3  =  1.933  +  jO. 697 

1.0  +  jO.O 

6 

-0.401  -  jO. 164 

a1  =  1.0  +  jO.O 
a2  =  -0.408  +  jO.  354 
a3  »  1.933  -  jO.  697 

0.000  +  jO.  051 
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Table  VII 


Stoneley  Wave  Solution  for  ±15°  Copper 
Velocity  =  1392.403  m/s 


N 

Roots 

bN 

Displacement  Vectors 

“l 

Weighting  Factors 

CN 

1 

0.0  -  jO. 115 

=  1.0  +  jO.O 

02  =  -76.18  +  jO.O 
o3  =  0.09  +  j  14.  011 

0.002  +  jO. 055 

-  -  -  --  --  - 

2 

0.0  +  jO. 115 

o^  =  1.0  +  jO.O 

02  =  76. 18  +  jO. 0 
o3  =  0.0  -  j  14.  011 

-0.002  -  jO. 055 

3 

+0.407  -  jO. 285 

o^  =  1.0  +  jO.O 
o2  -  -0.229  -  jO.  534 
o3  =  1.566  -  jO. 927 

-1.0  +  jO.O 

I 

+0.407  +  jO. 285 

-  ...... 

a3  =  1.0  +  jO.O 
a2  =  0.229  -  jO. 534 
a3  *  -1.566  +  jO. 927 

-0.998  +  jO. 055 

5 

-0.407  -  jO. 285 

*  1.0  +  jO.O 
o2  =  -0.  229  +  jO.  534 
o3  -  1.566  -  jO. 927 

0.998  -  jO. 055 

6 

-0.407  +  jO. 285 

oj  *  1.0  +  jO.O 
a2  *  0.  299  +  jO.  534 
a3  -  1.566  +  jO. 927 

1.0  +  jO.O 
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4.  Particle  motion  in  the  direction  of  wave  propagation  is  always 
zero  at  the  interface. 

These  characteristics  became  useful  guidelines  in  future  searches 
conducted  on  more  complex  crystal  systems.  Although  numerically  modified,  all 
except  the  last  were  found  to  be  true  in  all  cases  studied. 

4. 2  Non-Cubic  Materials 

Following  the  work  on  cubic  materials  described  in  Section  4. 1  a 
study  of  non-cubic  materials,  in  particular  quartz  and  LiNb03  was  begun. 

These  materials  are  of  particular  interest  because  they  are  the  more  commonly 
used  piezoelectric  substrates  for  surface  wave  devices. 

Our  first  attempt  was  to  search  along  major  axis  directions  where 
simple  surface  wave  solutions  were  well  known.  From  our  studies  of  simple 
cubic  materials  we  knew  that  Stoneley  wave  solutions  were  most  likely  to  be 
found  along  those  directions  where  the  normal  surface  acoustic  wave  existed 
with  only  t»*>  particle  motions,  both  in  the  saggital  plane.  In  the  case  of 
quartz  none  of  the  major  cuts,  X,  Y,  or  Z  yielded  Stoneley  wave  solutions. 

Failing  to  find  a  major  cut  of  quartz  which  would  support  a  Stoneley 
wave  our  attention  focused  upon  LiNbOj.  A  Stoneley  wave  search  was  carried 
out  on  the  three  major  cuts  of  this  material  and  only  the  Y-cut  was  found  to 
support  Stoneley  waves  for  propagation  along  the  Z  axis. 
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Stoneley  wave  solutions  on  YX  LiNbO-j  were  found  only  over  a  small 
range  of  material  angular  rotations,  i.e.,  a  maxim'im  of  +10  degrees.  Shown  in 
Table  VIII  is  a  typical  solution  for  ±5°.  Because  this  material  is  piezo¬ 
electric  eight  complex  decay  constants  are  shown  along  with  the  eight  partial 
solution  vectors.  The  total  displacements  and  potential  for  the  unshorted 
interface  are  plotted  in  Fig.  8.  Shown  in  Fig.  9  are  the  disp1  accents  and 
potential  when  the  interface  is  shorted  by  a  conducting  layer  of  zero  physical 
thickness.  In  this  case  the  velocity  change  was  comparable  to  an  electro¬ 
mechanical  coupling  of  1  *  10~^. 

An  important  characteristic  of  these  waves  which  was  substantially 
different  from  earlier  results  with  cubic  materials  was  the  low  decay  con¬ 
stants  This  is  shown  in  Fig.  9  where  clearly  more  than  25  wavelengths  are 
required  for  the  bound  solution  to  decay. 

At  the  present  time  our  studies  on  quartz  and  LiNbOj  are  incomplete, 
however,  Stoneley  waves  have  been  found  in  at  least  one  case.  The  effects  of 
piezoelectricity  upon  these  waves  are  not  well  understood  and  further  study  is 
needed.  Also  other  off-axis  orientations  may  yield  more  suitable  cuts  with 
larger  decay  constants  for  Stoneley  wave  devices. 
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Fig.  9  Shorted  interface  Stoneley  wave  displacements  and 
potential  for  ±5  ,  YZ  LiNbO^. 
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Table  VIII 

Stone! ey  Wave  Solution  for  ±2°  YZ  LiNb03 
Velocity  =  3527.576  m/s 
(Page  1  of  2) 


Roots 

bN 

Displacement  Vectors 

ai 

Weighting  Factors 

CN 

-0.  373  -  jl.  035 

aj  =  1.121  +  jl.  165 
a2  =  0.115  -  jO. 085 
a-j  =  1.0  +  jO.O 
a4  =  1.041  +  j  1.934 

-1.009  x  10'2  ♦  j2. 004  x  10’6 

-0.373  +  j2. 035 

aj  =  1.121  -  jl. 165 
a2  =  -0.115  -  jO. 085 

a-j  =  1.0  +  j  0.0 

a4  =  1.041  -  jl. 934 

-982  x  10‘3  -  j2. 315  x  10'3 

0.394  -  jO.  768 

=  0.664  +  jO. 651 
a2  =  0.160  +  jO.  043 
a-j  =  1.0  +  jO.O 
a4  =  0.673  +  jl. 592 

=2.061  x  10'2  -  j3. 31  *  10'4 

0.394  +  jO. 768 

=  -0.664  -  jO.  651 
a2  =  0. 160  +  jO. 043 
a3  =  1.0  +  jO.O 
a4  =  0.678  -  jl. 592 

-1.014  x  IQ'2  -  j4. 403  x  10'3 
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Table  VIII 

Stoneley  Wave  Solution  for  ±2°  YZ  LiNbC^ 
Velocity  =  3527.576  m/s 
(Page  2  of  2) 


Roots 

bN 

Displacement  Vectors 

ai 

Weighting  Factors 

CN 

0.106  -  jO. 134 

=  0.111  +  jO. 013 
a2  =  1.396  +  j  520 
a-j  =  1.0  +  jO.O 
a4  =  0.103  -  jO. 726 

-6.027  x  10'2  -  j 1. 264  x  lO*1 

0.106  +  jO. 134 

=  0.111  -  jO. 013 
a2  =  1.  396  +  j5. 52 
ct^  =  1.0  +  jO.O 
a4  =  0.103  +  jO. 726 

-8.764  x  10‘2  +  jl . 092  x  IQ'1 

-0. 057  -  jO. 009 

=  0.055  +  jO. 007 
a2  =  0.  554  +  j0.O25 
=  1.0  +  jO.O 
a4  =  0.010  -  jO. 0134 

9.734  x  10'1  +  j2. 292  x  IQ*1 

-0.057  +  jO. 009 

aj  =  0.055  -  jO. 007 
a2  =  0.  554  +  jO .  025 
=  1.0  +  jO.O 
a4  =  0.010  +  jO.  0134 

1.0  +  jO.  0 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  objective  of  this  study  was  to  assess  the  feasibility  of  develop' 
ing  a  single  crystal  pressure  sensor  utilizing  Stoneley  waves.  A  theoretical 
study  of  Stoneley  wave  propagation  was  carried  out  on  cubic  material  as  well 
as  non-cubic  but  piezoelectric  quartz  and  L i N bO 3 .  The  results  of  these 
studies  showed  that  strongly  bound  Stoneley  waves  could  propagate  along  major 
axes  of  cubic  materials.  Extending  the  analysis  to  the  non-cubic  materials, 
quartz  and  LiNbO^,  revealed  that  the  major  axes  of  quartz  could  not  support 
Stoneley  waves  while  only  the  YZ  cut  of  LiNbO-j  could  support  Stoneley  wave 
propagation.  The  analysis  for  LiNbOj  showed  that  the  waves  would  have  an 
electromechanical  coupling  coefficient  of  1  *  10*4  and  would  require  approx¬ 
imately  50  wavelengths  of  material  to  isolate  the  interface  from  the  outer 
surfaces  of  the  crystal  structure. 

In  view  of  the  level  of  effort  we  find  these  results  encouraging  but 
incomplete  in  the  sense  that  we  were  not  able  to  study  many  other  potential 
materials  and  orientations.  We  recommend  the  following  tasks  for  future  work 
in  this  area: 


1.  To  develop  a  better  understanding  of  the  role  played  by  crystal 
anisotropy,  surface  wave  behavior,  and  bulk  wave  properties  in 
Stoneley  wave  propagation. 

2.  To  utilize  this  knowledge  to  develop  better  and  less  time  con¬ 
suming  Stoneley  wave  search  methods. 

3.  To  examine  other  materials  such  as  silicon  and  GaAs  for  possible 
Stoneley  wave  propagation. 
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